Nonequilibrium resonant spectroscopy of molecular vibrons 
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Quantum transport through single molecules is essentially affected by molecular vibrations. We 
investigate the behavior of the prototype single-level model with intermediate electron- vibron cou- 
pling and arbitrary coupling to the leads. We have developed a theory which allows to explore 
this regime via the nonequilibrium Green function formalism. We show that the nonequilibrium 
resonant spectroscopy is able to determine the energies of molecular orbitals and the spectrum of 
molecular vibrations. Our results are relevant to scanning tunneling spectroscopy experiments, and 
demonstrate the importance of the systematic and self-consistent investigation of the effects of the 
vibronic dynamics onto the transport through single molecules. 
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Molecular electronics [H, 0, S| is the most promising 
development of nano-electronics, raising new theoretical 
challenges and calling for new methods. The understand- 
ing of quantum electron transport at the molecular scale 
is a key step to future devices. Recently, the interac- 
tion of electrons with molecular vibrations attracted an 
enormous attention unveiled by electron transport ex- 
periments through single molecules [1, 0, 0, 0, HI, 0] . In 
scanning tunneling spectroscopy (STS) experiments clear 
signatures of the electron-vibron interaction have been 
observed 0, ELI EH Il3j |. In these experiments both the 
electron-vibron coupling and the electron hybridization 
are well beyond the perturbation limit, so that a the- 
ory beyond the linear response or master equation ap- 
proaches is necessary. In this Letter, we present such a 
theory and discuss the vibronic features in the current- 
voltage curves of molecular junctions as a spectroscopic 
tool. 

Resonant electron transport through a molecular or- 
bital with electron-vibron interaction at finite voltages is 
essentially different from the well known inelastic electron 
tunneling spectroscopy (IETS). Indeed, in the IETS sit- 
uation one has typically the multi-channel tunneling be- 
tween two bulk metals with continuous density of states, 
and an increase of the conductance is observed when the 
bias voltage exceeds the threshold determined by the fre- 
quency of vibrational mode. As wc shall see below, in 
resonant tunneling through molecular levels coupled to 
vibrons, the situation is very different: in addition to the 
usual IETS signal, the inelastic signal can be observed 
as an extra peak in the differential conductance, or just 
as a decrease of the conductance. The resulting signal 
is determined by the energy of the electronic state, the 
vibronic frequency, the electron-vibron coupling, and the 
molecule-to-lead coupling. The position of the molecular 
energy levels could be changed by the external gate volt- 
age, but in the STS experiments the application of the 
gate voltage is not possible, instead the tip-to-molcculc 
distance can be changed, also changing the molecule-to- 
tip coupling llj . Thus it is very important to investigate 
quantum transport at arbitrary coupling to the leads. 
The simplest model, which describes properly these pe- 
culiarities, is a single electron with energy eo interacting 
with a single vibron with frequency loq (Fig. [I]) , and cou- 
pled to two noninteracting equilibrium electrodes. 



We analyze here the electron-vibron model developing 
the prelimi nary results obtained in our previous publi- 
cations [3, [Hf. We consider a molecule coupled to free 
conduction electrons in the leads by a usual tunneling 
Hamiltonian. Furthermore, the electrons are coupled to 
vibrational modes. The full Hamiltonian is the sum of 
the molecular Hamiltonian Hm , the Hamiltonians of the 
leads Hjt(Lj , the tunneling Hamiltonian Ht describing 
the molcculc-to-lead coupling, the vibron Hamiltonian 
Hy including electron-vibron interaction and coupling of 
vibrations to the environment (describing dissipation of 
vibrons) 



H = H 



M 



Hv + Ht. + Hn + Hi 



(1) 



A molecule is described by a set of localized states \a) 
with energies e a and inter-orbital overlap integrals t a p 
by the following model Hamiltonian: 



did 



0' 



(2) 



where d^ a ,d a are creation and annihilation operators in 
the states \a), and ip a {t) is the (self-consistent) electrical 
potential. The index a is used to mark single-electron 
states (atomic orbitals) including the spin degree of free- 
dom. The parameters of a tight-binding model could be 
determined by ah initio methods This is a compro- 




Figure 1: (Color online) Schematic picture of the considered 
fluctuating single-level model. 
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Figure 2: (Color online) Spectral function at different 
electron-vibron couplings: A/cjo = 0.4 (black), X/ujo = 1.2 
(blue, dashed), and A/u>o = 2 (red); at eo/wo = 5, Fl/uJo = 
Tij/wo = 0.1. In the insert the spectral function at A/wo = 1.2 
is shown at finite voltage, when the level is partially filled. 
Energies are in units of Hloq. 



mise, which allows us to consider complex molecules with 
a relatively simple model. 

The isolated single-level electron-vibron model is de- 
scribed by the Hamiltonian (h = 1) 

Hm+v = (eo + etpa)cftd + uj a a)a + A (at + a) d'd, (3) 

where the first and the second terms describe the free 
electron state and the free vibron, and the third term is 
electron-vibron minimal coupling interaction. 

The electrical potential of the molecule tpo plays an im- 
portant role in transport at finite voltages. It describes 
the shift of the molecular level by the bias voltage, which 
is divided between the left lead (tip), the right lead (sub- 
strate), and the molecule as ipo = ipn + t)(<pl — <Pr) 01- 

The Hamiltonians of the right (R) and left (L) leads 
read 



(4) 



<Pi(t) are the electrical potentials of the leads. Finally, 
the tunneling Hamiltonian 



£1 T = 



E 

i—L,R ka,a 



X] [Vika, a 4 ka d a +h.C. 



(5) 



describes the hopping between the leads and the 
molecule. A direct hopping between two leads is ne- 
glected. 

Though the model described above has a long his- 
tory, the many questions are not answered up to now. 
While the isolated electron-vibron model can be solved 
exactly by the so-called polaron or Lang-Firsov trans- 
formation, the coupling to the leads produces a true 
many-body problem. The inelastic resonant tunneling 
of electrons coupled to phonons was first considered in 



Refs. [17|, UM, |19|, [20|. There the exact solution in the 
single-particle approximation was derived, ignoring com- 
pletely the Fermi sea in the leads. At strong electron- 
vibron couplings and weak couplings to the leads the 
satellites of the main resonant peak are formed in the 
spectral function (Fig. [5]). The question which remains 
is whether these side-bands can be observed in the dif- 
ferential conductance. New theoretical treatments were 
presented in Refs. [HI, [H, 0, 0, 0, 
31, 32| ■ In parallel, the theory of inelastic resonant tun- 
neling in scanning tunneling spectroscopy was developed 



33, 34, 3 



The essential progress in calculation of transport prop- 
erties in the strong electron-vibron interaction limit has 
been made with the help of the master equation approach 
(23, 0, Hi]]. This method, however, is valid only in the 
limit of very weak moleculc-to-lead couplings and ne- 
glects all spectral effects, which are the most important 
at finite coupling to the leads. 

Our approach is based on the noncquilibrium Green 
function technique (3fj| . l37l . |38| , which is now a standard 
method in mesoscopic physics and molecular electronics. 
We follow the formulation pioneered by Meir, Wingreen, 
and Jauho [j^ 0, El ■ This method was also further 
developed in Refs. [2J 



27 



I2& 12a 132, 135(. The case of 
intermediate and strong electron-vibron couplings at fi- 
nite tunneling rates is the most interesting, but also the 
most difficult. Only the approaches by Flensberg (zj ], 
and Galperin, Ratner, and Nitzan [32^ exist, both start- 
ing from the exact solution for the isolated system and 
then switching on tunneling as a perturbation. Our ap- 
proach is exact for non-interacting systems with arbitrary 
tunneling coupling, and is approximate in the strong 
electron-vibron coupling limit without tunneling. Still, 
its advantage is clear in the truly nonequilibrium formu- 
lation. 

The current in the left (i = L) or right (i = R) contact 
to the molecule is described by the well-known expression 



Ji= 



i=L,R 



i - h j^{^-e, i )(G<(e) 
f^e-e Vi )[G R (e)-G A (e)})}, 



(6) 



where /f (e) is the equilibrium Fermi distribution function 
in the lead i with the chemical potential /i.; , and the level- 
width function r i=L(K) (e) = T. laf3 (e) is 



r M/3 (e) = 2-K^2v ika ^V* kaa 8{e - e lka ). 



(7) 



The lesser (retarded, advanced) Green function matrix 
of a nonequilibrium molecule G <( - R,A ^ = G^g can be 
found from the Dyson-Keldysh equations in the integral 
form 



G*(e) = G*(e) + G*(e)S*(e)G*(e), 
G<( £ ) = G^(e)S<( e )G^(e), 



(8) 
(9) 



or from the corresponding equations in the differential 
form (see Refs.ll4l Il5l and references therein). 
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Figure 3: Differential conductance of a symmetric junction 
(77 = 0.5, Tr — Tl) at different molecule-to-lead coupling, 
from Fl/^o — 0.1 (lower curve) to Fl/^o = 10 (upper curve), 
X/luq = 1, eo/uJo = 2. Voltage is in the units of fkJo- 
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Figure 4: Differential conductance of an asymmetric junction 
(77 = 0, Tr = 20Fl) at different molecule-to-lead coupling, 
from r_R/ix)o = 0.2 (lower curve) to F_r/cjo = 4 (upper curve), 
\/u>o — 2, eo/uoo = 5. The voltage is in the units of fiuja 



Here £ 



R,< _ 



-,-R,<(T) 



,fl,<(T) 
'R 



+ £«.<W is the 



total self-energy of the molecule composed of the tun- 
neling (coupling to the left and right leads) self-energies 



v fl,<(T) T R,<(T) 
Zj j=L,R ~ 



J2kcr { 



V, 



J l % Vjka,p}, and the 



_ V K,<(V) 



a/3 



Here we as- 



vibronic self-energy E- B > < ( vr J = £ 

sume, that vibrons are in equilibrium and are not excited 
by the current, so that the self-consistent Born approxi- 
mation is a good starting point. The corresponding self- 
energies can be found in Refs. 0, ES US EE US [H, US 

For the single-level model all equations are significantly 
simplified. Combining Ji and Jr the expression for 
the current can be written for energy independent ^l(r) 
(wide-band limit) as 



J 



H R 



271-/1 r 



R 



r. 



deA(e) [/°(e - e^ L ) - /°(e - e^)] . 

(10) 

It looks as simple as the Landaucr-Biittikcr formula, but 
it is not trivial, because the spectral density A(e) = 
—2lmG R (e) now depends on the distribution function of 
electrons inside the fluctuating molecule and hence the 
applied voltage, fL = —ipR — V/2 [3]. A typical exam- 
ple of the spectral function at zero voltage is shown in 
Fig. [2j At finite voltage it is calculated self-consistcntly. 

Before considering the experiments, let us discuss a 
general picture of the vibronic transport in symmet- 
ric and asymmetric molecular bridges coupled to the 
electrodes, provided in experiments by the MCBJ and 
STM technique, respectively. The differential conduc- 
tance, calculated at different molecule-to-lead coupling, 
is shown in Fig. [3] (symmetric) and Fig. 0] (asymmetric). 
At weak coupling, the vibronic side-band peaks are ob- 
served, reproducing the corresponding peaks in the spec- 
tral function (Fig. [2]). At strong couplings the broad- 
ening of the electronic state hides the side-bands, and 
new features become visible. In the symmetric junc- 



tion, a suppression of the conductance at V ~ ±hu)Q 
takes place as a result of inelastic scattering of the co- 
herently transformed from the left lead to the right lead 
electrons. In the asymmetric junction, the usual IETS 
increasing of the conductance is observed at a negative 
voltage V ~ —ftwQ, this feature is weak and can be ob- 
served only in the incoherent tail of the resonant con- 
ductance. We conclude, that the vibronic contribution 
to the conductance can be distinguished clearly in both 
coherent and tunneling limits. 

Now let us discuss the STS experiments US El E El- 
Here we concentrate mainly on the dependence on the 
tip-to-molecule distance [llj. When the tip (left lead in 
our notations) is far from the molecule, the junction is 
strongly asymmetric: <C T^j and 77 — s- 0, and the 
conductance is similar to that shown in Fig. [4] When the 
tip is close to the molecule, the junction is approximately 
symmetric: w Yr and 77 « 0.5, and the conductance 
curve is of the type shown in Fig. [3] We calculated the 
transformation of the conductance from the asymmetric 
to symmetric case (Fig. [5]). It is one new feature appeared 
in asymmetric case due to the fact that we started from 
a finite parameter 77 = 0.2 (in the Fig. [H77 = 0), namely 
a single peak at negative voltages, which is shifted to 
smaller voltage in the symmetric junction. The form and 
behavior of this peak is in agreement with experimental 
results [ll|. 

In conclusion, at weak molecule-to-lead (tip, substrate) 
coupling the usual vibronic side-band peaks in the dif- 
ferential conductance are observed; at stronger coupling 
to the leads (broadening) these peaks are transformed 
into step-like features. A vibronic-induced decreasing 
of the conductance with voltage is observed in high- 
conductance junctions. The usual IETS feature (increas- 
ing of the conductance) can be observed only in the case 
of low off-resonant conductance. By changing indepen- 
dently the bias voltage and the tip position, it is possible 
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Figure 5: (Color online) Differential conductance at different 
molecule-to-STM coupling (see the text), from asymmetric 
junction with Yl/^o = 0.025, Fr/ujo = 0.5 and r\ = 0.2 
(lower curve, blue thick line) to symmetric junction with 
Tl/loo = Tr/loo — 0.5 and 77 = 0.5 (upper curve, red thick 
line), X/ljo = 1, = 2. Voltage is in the units of hu>o 



to determine the energy of molecular orbitals and the 
spectrum of molecular vibrations. In the multi-level sys- 
tems with strong electron-electron interaction further ef- 
fects, such as Coulomb blockade and Kondo effect, could 
dominate over the physics which we address here; these 
effects have to be included in a subsequent step. 

We acknowledge fruitful discussions with Peter Hanggi. 
This work was supported by the Volkswagen Foundation 
under grant 1/78 340, by the EU under contract IST- 
2001-38951, and by the DFG CU 44/3-2. 



[1] M. A. Reed and J. M. Tour, Sci. Am. 282, 86 (2000). 
[2] C. Joachim, J. K. Gimzewski, and A. Aviram, Nature 

408, 541 (2000). 
[3] G. Cuniberti, G. Fagas, and K. Richter, Introducing 

Molecular Electronics (Springer- Verlag, 2005). 
[4] M. A. Reed, C. Zhou, C. J. Muller, T. P. Burgin, and 

J. M. Tour, Science 278, 252 (1997). 
[5] H. Park, J. Park, A. K. L. Lim, E. H. Anderson, A. P. 

Alivisatos, and P. L. McEuen, Nature 407, 57 (2000). 
[6] J. Park, A. N. Pasupathy, J. I. Goldsmit, C. Chang, 

Y. Yaish, J. R. Petta, M. Rinkoski, J. P. Sethna, H. D. 

Abrun, P. L. McEuen, et al, Nature 417, 722 (2002). 
[7] W. Liang, M. P. Shores, M. Bockrath, J. R. Long, and 

H. Park, Nature 417, 725 (2002). 
[8] R. H. M. Smit, Y. Noat, C. Untiedt, N. D. Lang, M. C. 

van Hemert, and J. M. van Ruitenbeek, Nature 419, 906 

(2002). 

[9] L. H. Yu and D. Natelson, Nano Lett 4, 79 (2004); L. H. 
Yu, Z. K. Keane, J. W. Ciszek, L. Cheng, M. P. Stewart, 
J. M. Tour, and D. Natelson, Phys. Rev. Lett. 93, 266802 
(2004). 

[10] X. H. Qiu, G. V. Nazin, and W. Ho, Phys. Rev. Lett. 92, 
206102 (2004). 

[11] S. W. Wu, G. V. Nazin, X. Chen, X. H. Qiu, and W. Ho, 



[12 
[13 
[14 
[15 
[16 

[17 

[18 

[19 

[20 
[21 

[22 

[23 

[24' 
[25' 

[26 

[27 

[28 

[29 

[30; 

[31 

[32 
[33 
[34 
[35 

[36 
[37 
[38 
[39 
[40 
[41 



Phys. Rev. Lett. 93, 236802 (2004). 

J. Repp, G. Meyer, S. M. Stojkovic, A. Gourdon, and 

C. Joachim, Phys. Rev. Lett. 94, 026803 (2005). 

J. Repp, G. Meyer, S. Paavilainen, F. E. Olsson, and 
M. Persson, Phys. Rev. Lett. 95, 225503 (2005). 

D. A. Ryndyk and J. Keller, Phys. Rev. B 71, 073305 
(2005). 

D. A. Ryndyk, M. Hartung, and G. Cuniberti, Phys. Rev. 
B 73, 045420 (2006). 

S. Datta, W. Tian, S. Hong, R. Reifenberger, J. I. Hen- 
derson, and C. P. Kubiak, Phys. Rev. Lett. 79, 2530 
(1997). 

L. I. Glazman and R. I. Shekhter, Zh. Eksp. Teor. Fiz. 
94, 292 (1988), [Sov. Phys. JETP 67, 163 (1988)]. 
N. S. Wingreen, K. W. Jacobsen, and J. W. Wilkins, 
Phys. Rev. Lett. 61, 1396 (1988). 

N. S. Wingreen, K. W. Jacobsen, and J. W. Wilkins, 

Phys. Rev. B 40, 11834 (1989). 

M. Jonson, Phys. Rev. B 39, 5924 (1989). 

U. Lundin and R. H. McKenzie, Phys. Rev. B 66, 075303 

(2002) . 

J.-X. Zhu and A. V. Balatsky, Phys. Rev. B 67, 165326 

(2003) . 

S. Braig and K. Flensberg, Phys. Rev. B 68, 205324 
(2003). 

K. Flensberg, Phys. Rev. B 68, 205323 (2003). 

V. Aji, J. E. Moore, and C. M. Varma, cond-mat/0302222 

(2003). 

A. Mitra, I. Aleiner, and A. J. Millis, Phys. Rev. B 69, 
245302 (2004). 

T. Frederiksen, Master's thesis, Technical University of 
Denmark (2004). 

T. Frederiksen, M. Brandbyge, N. Lorente, and A.-P. 
Jauho, Phys. Rev. Lett. 93, 256601 (2004). 
M. Galperin, M. A. Ratner, and A. Nitzan, Nano Lett 4, 
1605 (2004); J. Phys. Chem. B 121, 11965 (2004); Nano 
Lett 5, 125 (2005). 

M. Cfzek, M. Thoss, and W. Domcke, Phys. Rev. B 70, 
125406 (2005). 

J. Koch and F. von Oppen, Phys. Rev. Lett. 94, 206804 
(2005); J. Koch, M. Semmelhack, F. von Oppen, and A. 
Nitzan, Phys. Rev. B 73, 155306 (2006). 
M. Galperin, M. A. Ratner, and A. Nitzan, Phys. Rev. 
B 73, 045314 (2006). 

B. N. J. Persson and A. Baratoff, Phys. Rev. Lett. 59, 
339 (1987). 

M. A. Gata and P. R. Antoniewicz, Phys. Rev. B 47, 
13797 (1993). 

T. Mii, S. G. Tikhodeev, and H. Ueba, Surf. Sci. 493, 
63 (2001); Surf. Sci. 502, 26 (2002); Phys. Rev. B 68, 
205406 (2003); S. G. Tikhodeev, and H. Ueba, Phys. Rev. 
B 70, 125414 (2004). 

L. Kadanoff and G. Baym, Quantum Statistical Mechan- 
ics (Benjamin, New York, 1962). 

L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964), [Sov. 
Phys. JETP 20, 1018 (1965)]. 

J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 
(1986). 

Y. Meir and N. S. Wingreen, Phys. Rev. Lett. 68, 2512 
(1992). 

A.-P. Jauho, N. S. Wingreen, and Y. Meir, Phys. Rev. B 
50, 5528 (1994). 

H. Haug and A.-P. Jauho, Quantum Kinetics and Optics 
of Semiconductors, vol. 123 of Springer Series in Solid- 
State Sciences (Springer, 1996). 



